THIS  REpORT  HAS  BEEK  DELIMITED 
AMD  CltiAREO  FOR  PUBLIC  RELEASE 

UNDER  DOD  DIRECTIVE  5200,20  AND 

NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  RELEASE,* 
DICTRIBUTION  UNLIMITED. 


ARP  A  Order  No.  529 


OJ 

CO 

00 


SURFACE  WAVES  ON  A  UNIAXIAL  PLASMA  SLAE 
THEIR  GROUP  VELOCITY  AND  POWER  FLOW 

by 

H.  L.  Bertoni  and  A.  Hessel 

Research  Report  No.  PIBMRI-1294-65 
Contract  No.  NONR  839(38) 

Sponsored  by 

Advanced  Research  Projects  Agency 

Monitored  by 
Office  of  Naval  Research 
Washington,  D.  C. 


October  12. 1965 


POLYTECHNIC  INSTITUTE  OF  BROOKLYN 
MICROWAVE  RESEARCH  INSTITUTE 


Electrophysic c  Department 


ARPA  Order  No.  529 


Research  Report  No.  PIBMRI- 1  294-65 
Contract  No.  NONR  839(38) 


SURFACE  WAVES  ON  A  UNIAXIAL  PLASMA  SLAB; 
THEIR  GROUP  VELOCITY  AND  POWER  FLOW 

by 

H.  L.  Bertoni  and  A.  Hessel 


Polytechnic  Institute  of  Brooklyn 
Electrophysics  Department 

Graduate  Center 
Route  110 

Farmingdale,  New  York  11735 


Research  Report  No.  PIBMRI-1294-65 
Contract  No.  NONR  839(38) 

October  12,  1965 


Title  Page 

Acknowledgement 

Summary 

Table  of  Contents 

20  Pages  of  Text 

References 


H.  L.  Bertoni 


A.  Hessel  A  ~ 

'\v4jfc  \ 

Herbert  J.  Carlin 
Head  of  Department 

A 


Reproduction  in  whole  or  in  part  is  permitted  for 
any  purpose  of  the  United  States  Government. 


Sponsored  by 

Advanced  Research  Projects  Agency 

Monitored  by 

Office  of  Naval  Research 
Washington,  D.  C. 


ACKNOWLEDGEMENT 


Thi  *  research  was  supported  by  the  Office  of  Naval  Research,  Washington,  D  C 

under  Contract  No.  NONR  839(38).  ARPA  Order  No.  529, 


SUMMARY 


The  relation  between  group  velocity  and  the  velocity  of  energy  transport  for  sur¬ 
face  waves  in  plane- stratified,  anisotropic,  dispersive  media,  which  was  derived  in  (1), 
is  verified  by  direct  calculation  for  the  case  of  surface  waves  on  a  uniaxial,  cold  plasma 
slab  located  in  free  space.  A  superimposed  D.  C.  magnetic  field  of  infinite  strength  and 
parallel  to  the  interfaces  generates  the  uniaxial  anisotropy  in  the  slab.  Surface  waves 
having  an  arbitrary  direction  of  propagation  with  respect  to  the  D.  C.  magnetic  field  are 
considered.  A  useful  graphical  presentation  of  the  dispersion  relation  is  given,  from  which 
the  direction  of  propagation  of  "surface  wave  rays"  is  directly  obtained. 
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INTRODUCTION 

The  group  velocity  of  a  surface  wave  propagating  in  a  plane- stratified,  anisotropic, 
dispersive  medium  that  is  also  linear  and  lossless  has  been  shown  to  be  equal  to  the  veloc¬ 
ity  of  energy  transport  of  the  surface  wave  as  a  whole.  ^  This  velocity  is  defined  as  the 
integral  of  the  real  part  of  the  Poynting  vector  over  the  coordinate  in  the  direction  of 
stratification  divided  by  the  integral  of  the  stored  energy  density  over  this  coordinate.  In 
this  paper  the  atave  relation  is  verified  by  direct  calculation  for  the  case  of  surface  waves 
supported  by  a  uniaxial,  cold-electron  plasma  slab.  The  plasma  slab  is  assumed  to  be  of 
infinite  extent  and  to  be  located  in  free  space.  A  static  magnetic  field  of  infinite  strength 
and  parallel  to  the  interfaces  between  the  plasma  and  free  space  generates  the  anisotropy. 

The  characteristics  of  trapped  surface  waves  propagating  on  anisotropic  plasma 

slabs  have  been  discussed  in  the  literature  for  various  specific  directions  of  propagation 

(2) 

relative  to  the  static  magnetic  field.  Wait  has  considered  the  surface  waves  propagating 

on  a  thin  plasma  slab  with  an  arbitrary  static  magnetic  field.  Requiring  the  plasma  slab 

to  be  thin  reduces  the  effect  of  the  static  magnetic  field  to  that  which  would  be  produced  by 

(3) 

the  component  normal  to  the  slab  alone.  Meltz  and  Shore'  discuss  the  excitation  of  sur¬ 
face  waves  on  a  slab  of  arbitrary  thickness  when  the  static  magnetic  field  is  perpendicular 
to  the  slab  and  of  infinite  strength.  In  both  of  these  cases,  the  anisotropy  is  such  that  the 
slab  configurations  have  rotational  symmetry  about  the  coordinate  normal  to  the  slab  and 
hence  the  characteristics  of  the  surface  waves  will  be  independent  of  the  direction  of  propa¬ 
gation.  Furthermore,  as  in  isotropic  slab  configurations,  the  velocity  of  energy  transport 
of  surface  waves  on  these  slab  configurations  will  be  parallel  to  the  transverse  wave  vector. 

When  the  static  magnetic  field  is  parallel  to  the  air-plasma  interfaces,  the  effect 
of  the  resultant  anisotropy  is  more  striking  since  then  the  characteristics  of  the  surface 
waves  on  the  slab  depend  on  their  directions  of  propagation  with  respect  to  the  static  mag¬ 
netic  field.  Also,  the  velocity  of  energy  transport  will  not  be  parallel,  in  general,  to  the 
transverse  wave  vector.  Examples  found  in  the  literature,  of  surface  waves  on  slab  con¬ 
figurations  with  axis  of  anisotropy  parallel  to  the  interfaces,  do  not  illustrate  these  aniso¬ 
tropic  effects  as  they  are  restricted  either  to  propagation  along^’  ^  or  normal^  to  the 
static  magnetic  field.  In  either  case,  the  velocity  of  energy  transport  is  parallel  to  the 
transverse  wave  vector.  In  this  paper,  however,  the  surface  wave  fields  are  considered 
for  arbitrary  directions  of  propagation  with  respect  to  the  static  magnetic  field  of  infinite 
strength.  It  will  be  shown  that  for  this  configuration,  the  velocity  of  energy  transport  of 
each  surface  wave  is  not  parallel,  in  general,  to  the  transverse  wave  vector,  and  that  the 
direction,  as  well  as  the  magnitude,  of  the  real  part  of  the  complex  Poynting  vector  varies 


2 


with  the  coordinate  normal  to  the  slab.  Thus  the  slab  configuration  provides  a  non-trivial 
example  of  the  equality  of  the  surface  wave's  group  velocity  and  its  energy  transport 
velocity.  The  excitation  of  the  surface  waves  is  act  considered  here. 

In  the  first  section  of  this  paper  the  fields  and  dispersion  relation  of  the  E-type 
surface  waves,  which  have  no  component  of  R.  F.  magnetic  field  along  the  static  magnetic 
field,  are  found.  A  graphical  procedure  for  solving  the  dispersion  relation  and  the  prop¬ 
erties  of  the  dispersion  curves  are  discussed  in  the  second  section.  The  last  section  is 
devoted  to  an  analytical  verification  of  the  equality  of  group  velocity  and  energy  transport 
velocity  for  the  surface  waves.  The  Appendix  contains  a  proof  that  the  uniaxial  slab  con¬ 
figuration  can  support  only  the  E-type  surface  waves  described  in  the  body  of  the  paper. 


FIELDS  AND  DISPERSION  RELATION 

In  this  section,  the  fields  and  dispersion  relation  for  surface  waves  on  a  uniaxial 
electron  plasma  slab  are  found.  The  plasma  within  the  slab  is  homogeneous  and  the  super¬ 
imposed  D.  C.  magnetic  field,  which  is  assumed  to  be  of  infinite  strength,  is  parallel  to 
the  y  axis  (see  Fig.  1).  In  the  linear  or  small  signed  approximation,  the  interaction  of 
the  uniaxial  plasma  with  a  monochromatic  electromagnetic  field  may  be  described  by  a 
relative  dielectric  tensor  c  '  .  Neglecting  collision  loss,  when  the  D.  C.  magnetic  field 
is  in  the  y  direction,  e'  takes  the  form  '  ^ 

’  1  0  o- 

C  e*  3  =  0  I-X  0  (1) 

-0  0  1. 


where 
slab  e 


2 

X  =  (uu  /it)'  and  u  is  the  electron  plasma  frequency.  Thus  in  the  plasma 
P  P 

=  eQ  e  '  while  in  the  air  regions  e  =  e.  1  ,  where  1  is  the  unit  dyadic.  The  per¬ 


meability  tensor  ia  is  given  everywhere  by  y  =  1  . 


The  surface  wave  fields,  which  decay  exponentially  in  the  air  regions,  have  trans- 
-j(kxx  +  k  y) 

verse  dependence  e  ’  ,  k^  and  k^  being  real  transverse  wave  numbers.  These 

fields  will  be  constructed  from  those  plane  wave  solutions  appropriate  to  the  plasma  reg¬ 


ion  and  those  appropriate  to  the  free  space  regions.  The  plane  wave  solutions  appropriate 


to  the  plasma  slab  are  those  waves  of  the  form 


E(r,k)| 

e'(k) 

H(r,k)l 

h'(k) 

(2) 


Fig.  1  Physical  Configuration  of  the  Anisotropic  Plasma  Slab. 


which  can  propagate  in  an  infinite  homogeneous  plasma  described  by  the  relative  dielectric 
tensor  e'  given  in  (l).  Similarly,  the  plane  wave  solutions  appropriate  to  the  free  space 
are  those  having  the  form  given  in  (2)  which  can  exist  when  the  plasma  slab  is  absent. 

The  plasma  plane  waves  are  found  by  substituting  IS  and  H  from  (2)  into  Maxwell's 
equations.  The  resultant  equations  are,  when  the  time  dependence  e^  ^  ;s  suppressed, 


kite'  =  jj  u  Q  h  ' 
k  xh'  =  -  ju  cQ  e'-  e  ' 


(3) 


2 

Multiplying  the  first  equation  by  kx  and  substituting  the  second  gives  kx(k  x  e')=-k  «'• 

2  2  ,  . 
with  kQ  =  uu  eoUQ.  Expanding  the  triple  cross-product,  this  equation  maybe  written  in 

dyadic  form  as 


e, 


(k^c'+ kk-k2  1)  •  e'  =0 


(4) 


which  is  equivalent  to  three  homogeneous  equations  in  three  unknowns.  For  there  to  be 

non-trivial  solutions  of  (4),  the  determinant  of  the  matrix  representation  of  the  dyadic 
2  2 

operator  (kQc'+kk-k  1)  must  vanish.  Letting  k=kt+_zoH>  with  k{  =  x ok^  +  ^0^y’  t^le 
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vanishing  of  the  determinant  results  in  the  plane  wave  dispersion  relation  D^(k ,  (j)  =  0 
for  an  infinitely  extended,  homogeneous,  uniaxial  plasma.  This  plane  wave  dispersion 
relation  maybe  solved  for  h-  Four  solutions  result,  which  are 


The  sign  choice  before  each  root  refers  to  waves  carrying  power  or  decaying  in  the 
positive  or  negative  z  direction.  Substituting  each  of  the  four  solutions  given  in  (5)  into 
(4),  the  corresponding  field  vector  e'  can  be  found.  Finally,  the  pertinent  h'  can  be 
calculated  from  (3). 

The  above  method  may  be  repeated  to  find  the  plane  wave  fields  e  '  and  h'  for 
free  space.  Since  kt  must  be  the  same  for  the  entire  urface  wave  if  the  transverse 
fields  are  to  be  continuous  everywhere  across  the  planes  z  =  ±  d,  it  follows  that  the  free 
space  wave  vector  k  =  k  +z  h  .  Substituting  the  form  of  E  and  H  given  in  (2)  into 
Maxwell's  equations  for  free  space  gives 


k  x  e  '  =  Oiii  h' 
—  a  —  a  o — a 


k  xh1  -  -uu  e  e  ' 
—  a  —  a  o—  a 


(6) 


From  (6), 
be 


the  homogeneous  equations  that  determine 


e 

—  a 


are  found,  in  dyadic  form,  to 


[<k2-k2)  1  +  k  k  ]•  e' 
o  —a  ~  — a— aJ  —a 


(7) 


From  the  requirement  that  the  determinant  of  the  matrix  representation  of 
2  2 

[(k  -k  )  1+  k  k  ]  vanish  for  non-trivial  solutions  of  e'  to  exist,  one  can  solve  for 
L  o  —  a  ~  — a— aJ  —a 

xa  as 

When  these  values  of  k  are  used,  (7)  reduces  to  k  •  e  '  =  0,  i.  e.  ,  the  plane  wave 

a  —  a  —  a 

electric  field  is  orthogonal  to  the  wave  vector  k  ,  a  condition  that  does  not  uniquely 

a 

determine  e#  .  Commonly  chosen  solutions  for  e/  are  those  corresponding  to  TM  and 
*“  a  a 

TE  modes  with  respect  to  the  z  direction.  Other  possible  choices  for  e^f  which  will 
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prove  more  useful  in  this  analysis,  are  those  of  the  so-called  E-type  and  H-type  modes, 
which  are  appropriate  linear  combinations  of  the  TM  and  TE  modes.  The  E-type  modes 
with  respect  to  y  are  characterized  by  the  vanishing  of  the  y  component  of  the  magnetic 
field,  while  the  H-type  modes  are  characterized  by  the  vanishing  of  the  y  component  of 
the  electric  field. 

In  each  region,  the  surface  wave  fields  will  be  a  combination  of  the  plane  wave 
solutions  appropriate  to  that  region,  the  relative  amplitudes  of  which  can  be  found  from 
the  radiation  condition  and  the  continuity  conditions  at  z  =  ±d.  Since  the  free  space 
outside  the  slab  is  homogeneous,  the  surface  waves  are  characterized  by  an  exponential 
decay  of  their  fields  away  from  the  slab.  Such  decay  requires  that  h  be  imaginary  and 

Si 

that  for  Region  1  the  sign  choice  in  (8)  be  taken  to  give  n  =  -j|n  |  ,  so  that  the  fields  will 

decay  in  the  positive  z  direction.  For  the  plane  waves  in  Region  3,  the  sign  must  be 

taken  so  as  to  give  k  =  j |  k  | ,  which  will  result  in  fields  that  decay  in  the  negative  z 

3i  Si 

direction. 


It  will  be  shown  later  that  a  surface  wave,  whose  fields  in  the  plasma  slab  are  a 

i  — 2  2  T~ 

combination  of  those  plane  wave  fields  corresponding  to  h=  ±V(  1 -X)(kQ  -  ky)  -  k^  ,  exists 
only  for  X>1.  In  the  Appendix  it  is  shown  that  H-type  surface  wave  modes,  characterized 
by  the  vanishing  of  the  y  component  of  the  electric  field,  cannot  propagate  on  the  uniaxial 

plasma  slab.  The  plane  wave  fields  in  the  plasma  region  corresponding  to 

r - ~7 — 2  ~'l~ 

x=  ±V(1-X)(k  -k  )-k  are  E-type  modes  and  have  the  form 


e '  =  A[x  k  k 
—  —  o  x  y 


y  (k^  -  k^)  +  z  k  h  ] 
-t-o  o  y  — o  y 


h'  =  A  W€o[xoH-  zokxl 
with  A  an  arbitrary  constant. 

A s  mentioned  earlier,  the  only  requirement  on  e  '  is  that  k 


—a 


—a  —a 


(9) 


e  '  =  0.  Hence, 


we  may  arbitrarily  select  the  transverse  part  of  e^  and  then  use  the  requirement 
-a'  -a  =  °  t0  find  the  corresPonding  z  component  of  e^.  A  particularly  useful  form  of 
the  transverse  part  of  is  obtained  by  choosing  it  to  be  identical  with  the  transverse 
part  of  e'  as  given  in  (9).  This  choice  will  be  seen  to  simplify  the  application  of  the  con¬ 
tinuity  requirement  on  the  transverse  fields  at  z-  ±d.  Following  this  procedure  one  finds 

e  '  =  B^S0kxky  -  lJko"  ky)  +  -okyHa^ 

—  a  B  :t)  eo^— o  Ka  "  — o  kx^ 


(10) 
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with  B  an  arbitrary  constant.  It  is  seen  that  the  transverse  part  of  h^  has  the  same 
vector  direction  as  the  transverse  part  of  h It  will  thus  be  possible  to  satisfy  the  con¬ 
tinuity  conditions  at  z  -  id  us'ng  only  the  plane  waves  exhibited  in  (9)  and  (10). 


Since  the  w  wc  number  v  must  be  imaginary,  let 

a 

a  -  Vk2-k2 


(ID 


so  that  h  -  ±ja  where  a  is  real  and  positive.  Then  in  Region  1,  x  =  -jr:,  for 
a  0  a  * 

decay  in  the  positive  z  direction,  and  if  p=xQx+yoy,  the  fields  are 

£  s  Bl[iokxky'rJko'ky|--‘.ojkyal*’a*'  ^  ~ 

-jk.p  021 

H  "  B  i  -i  e  [-x  j  a  -  z  k  ]  e  3  z  e  - 

while  in  Region  3,  y  =  ja*  and  the  corresponding  fields  are 
a 

E=  B3^0kxky  -Io(k^-k^)  +_z<jjkya]eaZe  J_t  " 

-jk-o  •  (13) 

H  =  Biiefx  ja-z  k  ]  eaz  e  - 

—  J  O  — O  — O  X 


The  constants  and  have  yet  to  be  determined. 

For  simplicity  in  what  follows,  we  define  3  as 

3  =  V(  1-X)(k2-  k2  )  -  k2  (14) 

’  o  y  x  '  ' 

so  that  in  (5)  h  =  ±3-  It  will  be  shown  that  for  a  surface  wave  to  propagate  on  the  slab  P 
must  be  real.  The  fields  in  Region  2  will  be  the  sum  of  the  fields  of  the  two  plane  waves 
having  the  vector  form  displayed  in  (9)  and  traveling  in  opposite  directions  along  z.  The 
most  genera]  form  of  such  a  sum  is 


E  =  j[x  k  k  -  y  (k2-k2)](A  e  z  +  A,e^  8  z) 
—  I  — o  x  y  -k-  o'  o  y  J  1  2 


+  z  k  3  ( A .  e  °  Z  -  A^e^  Z  )  l  e  t- 
-o  v  1  2  J 


( 1 5  -  a) 


H  =  ■Juec{xo9(A1  e’j8z- A,e->8z) 


-z  k  (A  e"->3z  +  A,ej9z)le  J~t  - 
-oxl  2  J 


( 15-b) 
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t-:.- 


with  A^  and  A^  to  be  determined  from  the  boundary  conditions  at  z  -  ±d. 

Requiring  and  Ht  to  be  continuous  at  z  =  ±d  results  in  four  homogeneous 
equations  in  four  unknowns  from  which  the  relative  amplitudes  as  well  as  the  surface 
wave  dispersion  relation  can  be  found.  The  continuity  conditions  at  z  =  d  given  the  equa¬ 
tions 


B :  e”ad  =  Aie^Sd  +  A2ej8d 


"j  aB1  e*ad  =  B(Aie-j0d-A2ej0d) 


(16-a) 


while  those  at  z  =  -d  result  in 


B3e'ad  =  Aiej8d+  A2e^0d 
jaB3e‘ad=  0(A1eif3d-A2e-iBd) 


(16-b) 


Elimination  of  B, 
set 

0  = 

0  = 


from  the  first  two  equations  and  B, 

A  e'jM(l  +  -£-)  +  A?ej8d(l  -  J-) 

1  ja  2  ja 

A  J0d(l~2-)  +A-e‘j0d(i  +  ± ) 

1  ja  2  j  a ' 


f  .'om  the  second  two  gives  the 


(17) 


which  has  a  non- trivial  solution  for  A^  and  A.,  only  if  the  determinant  of  the  coefficients 
is  zero.  The  vanishing  of  the  determinant  yields  the  surface  wave  dispersion  relation 


»i48d-(j^f)2=°  •  <‘8> 

If  expressions  (11)  and  (14)  for  a  and  8  in  terms  of  k^,  k^  and  kQ  =  are  sub¬ 

stituted  into  (18),  the  dispersion  relation  is  seen  to  be  of  the  form  D  (k  ,  id)  =  0. 

S  t 

PROPERTIES  OF  THE  DISPERSION  RELATION 

As  given  in  (14),  8  is  either  real  or  imaginary  for  all  real  k^  and  k^.  Let  us 
first  verify  that  no  solutions  of  (18)  exist  for  which  8  is  imaginary.  If  8  is  imaginary, 
i.  e. ,  8  =  ±  j  1 8  1  >  then  (18)  becomes 


8 


e  +  4  |B  |d  /oililV  _  |19) 

The  left-hand  side  is  less  (greater)  than  unity  while  the  right-hand  side  is  greater  (less) 

than  unity.  This  contradiction  verifies  the  assertion.  When  0  is  real,  however,  both 

terms  in  (18)  have  magnitude  unity  so  that  a  solution  is  possible.  In  order  to  find  the 

range  of  frequencies  for  which  (18)  has  solutions,  based  on  the  restriction  that  0  be  real, 

we  plot  for  all  X  >  0  those  regions  in  the  k  -  k  plane  where  0  is  real  and  the  region 

x  y 

where  a  is  real  (see  Fig.  2).  From  Fig.  2  it  is  .,een  that  the  regions  where  0  is  real 
and  the  region  where  a  is  real  overlap  only  when  X  >1.  Hence  the  possibility  that  sur¬ 
face  waves  can  propagate  exists  only  for  X  >  1.  In  passing,  observe  that  (18)  remains 
invariant  under  the  substitution  of  -8  for  0  .  Thus  it  is  sufficient  to  consider  only  posi¬ 
tive  values  of  0.  Since  the  slab  configuration  has  mirror  symmetry  in  the  plane  z  =  0, 
the  surface  wave  fields  will  correspond  to  either  an  open-circuit  or  a  short-circuit 
bisection  of  the  slab  (even  and  odd  solutions  in  z).  The  dispersion  relation  given  in  (18) 
can  be  split  into  two  independent  dispersion  relations,  one  giving  the  open-circuit  bisection 
solutions  and  the  other  the  short-circuit  bisection  solutions.  These  are 


_j  2  0d  ..  .  j  a  +  8 
e  -  ‘jTCs 


(20) 


where  the  plus  and  minus  signs  correspond  to  short-circuit  and  open-circuit  bisections, 
respectively.  Using  the  plus  sign  for  the  short-circuit  bisection  case,  the  dispersion 
relation  may  be  put  in  the  form 

a  =  -8  cot  0d  ,  (21) 


whereas  if  the  minus  sign  is  used,  the  dispersion  relation  for  the  open-circuit  bisection 
case  can  be  written 

a  =  0  tan  0d  (22) 


with  a  and  0  as  given  in  (11)  and  (14). 

A  graphical  method  for  solving  equations  (21)  and  (22)  is  described  below. 

In  order  to  show  that  equations  (21)  and  (22)  are  satisfied  for  real  values  of  k  ,  k  and 

x  y 

tu  <  (Dp  .  i.  e.  ,  X  >  1,  the  plots  of  (21)  and  (22)  in  the  0  -  a  plane  are  considered. 
Since  a  and  0  have  been  taken  to  be  positive,  only  the  first  quadrant  is  of  interest. 


a  REAL  OUTSIDE 
CIRCLE  kj+lykj 


p-  0  FOR  X  >1 


PmO  FOR  X<l 


£*0  FOR  X>l 


Fig.  2  Regions  of  Real  S  and  a  in  the  kx-  ky  Plane. 
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2  2  2  2  2  2 

Adding  a  to  0  gives  a  +  0  =  X(k  -k  )  when  (11)  and  (14)  are  used.  Because 

2  2  y  0 

ky  >  kQ  for  surface  waves,  as  can  be  seen  from  Fig.  2,  the  plot  of  this  relation  in  the 
6  -  a  plane  is  a  circle  whose  radius  is  Vx o?  -  k^)  .  The  intersection  of  this  circle  with 
the  plot  of  (21)  or  (22)  gives  a  and  0  from  which  k^  can  be  found  using 

kx  =  UX-l)a2-02]  •  (23) 

But  kx  must  be  real  so  that  in  the  first  quadrant  only  those  intersections  for  which 
a  >-===■  0  give  values  of  a  and  3  which  correspond  to  an  actual  surface  wave.  Since 

2  2 

a  and  0  depend  only  on  k  andk  and  not  merely  on  k  ork  ,  constant  uu  surface  wave  dis- 

y  x  ’  y  x 

persion  curves  will  have  mirror  symmetry  about  the  k  and  k  axes  in  the  k  -k  plane. 

x  y  x  y 

Thus,  knowing  the  relation  between  k  andk  for  k  ,k  >0  is  sufficient  to  determine 

x  y  x  y 

the  entire  dispersion  curve. 

Figure  3  has  been  sketched  to  show  the  method  outlined  above  for  finding  that  k 

which  satisfies  (21)  when  k^  and  kQ  are  given.  Each  branch  of  -0cot0d  depicted  in 

Fig.  3  corresponds  to  a  particular  short-circuit  bisection  surface  wave  mode.  Since  there 

are  an  infinite  number  of  such  branches,  there  will  be  an  infinite  number  of  short-circuit 

2  2  2  2 

bisection  surface  wave  modes  when  X(k  -k  )  -.  ® .  For  a  finite  value  of  X(k  -k  )  only 

y  o  y  o  7 

a  finite  number  of  surface  wave  modes  can  propagate.  For  values  of  a  and  0  in  the 
shaded  region  of  Fig.  3,  k  ,  as  found  from  k  :  iJiL(X-l)a  -0  ],  is  imaginary.  Thu" 

X  XiA  ^ 

it  is  seen  that  for  fixed  k  each  mode  has  a  minimum  value  of  k^  >  k  at  which  k  =  0 

o  y  o  2  x 

and  below  which  no  real  solutions  for  k  exist.  The  minimum  value  of  k  for  which  a 

x  y 

particular  surface  wave  mode  can  exist  is  found  from  the  condition  that  the  circle 
2  2  2  2  1 

a  +  0  =  X(k  -k  ),  the  line  a  =  — - 0  and  that  branch  of  a  =  -0  cot  0d  corresponding 


to  the  mode  in  question  all  intersect  at  a  common  point.  As  k^  increases  from  its 
minimum  value,  k^  and  the  corresponding  solutions  for  k^  for  each  branch  of  -0cot0d 
trace  out  the  surface  wave  dispersion  curves  in  the  k  -  k  plane  of  the  short-circuit 
bisection  modes. 

In  a  similar  fashion,  Fig.  4  depicts  the  method  for  finding  that  k^  which  satisfies 
(22)  when  k^  and  kQ  are  given.  From  this  figure  and  Fig.  3,  it  is  seen  that  the  lowest 
surface  wave  mode  on  the  slab,  i.  e.  ,  the  one  with  the  smallest  value  of  0,  is  that  open- 
circuit  bisection  mode  corresponding  to  the  branch  of  0  tern  0d  starting  at  0=0.  As  in 
the  case  of  the  short-circuit  bisection  modes,  corresponding  to  values  of  0  and  a 
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in  the  shaded  region  of  Fig.  4  is  imaginary.  Thus  for  each  of  the  higher  open-circuit 

2  2 

bisection  modes  there  will  be  a  minimum  value  of  k  >k  at  which  k  =0  and  below 

y  o  x 

which  no  real  solution  for  k  exists.  For  the  lowest  open-circuit  bisection  mode 
d  ^ 

[^■(0  tan  0d)]g_Q  =  0  so  that  a  part  of  the  branch  of  0tan0d  starting  at  0  =  0  lies  in  the 

shaded  region  of  Fig.  4.  Hence  there  will  also  be  a  minimum  value  of  k^  >k^  for  the 

lowest  surface  wave  mode  below  which  no  real  solution  for  k  exists.  As  in  the  short- 

2  .  x 

circuit  bisection  case,  when  k^  increases  from  its  minimum  value  for  a  particular  mode 
and  for  a  fixed  kQ,  k  and  the  corresponding  value  of  k  trace  out  the  dispersion  curve 
of  that  open-circuit  bisection  mode. 

In  what  follows,  the  basic  properties  of  the  surface  wave  dispersion  curves  will 
be  derived.  For  any  one  mode,  these  properties  lead  to  the  form  of  the  dispersion 
curves  shown  in  Fig.  5,  which  has  been  drawn  for  two  different  frequencies  uu^  >  OJ^.  In 
order  to  find  the  shape  of  the  dispersion  curves  of  any  one  mode  and  for  fixed  uu,  we  con¬ 
sider  the  corresponding  branch  of  -0  cot  0d  in  Fig.  3  or  of  0  tan  0d  in  Fig.  4.  As 
pointed  out  previously,  the  dispersion  curves  are  symmetric  about  the  k^  and  k^  axes 
so  that  we  need  find  only  that  portion  of  the  curves  in  the  first  quadrant  of  Fig.  5.  Also, 
as  was  previously  discussed,  in  the  first  quadrant  of  Fig.  5,  k^  takes  on  its  minimum 
value,  which  is  greater  than  kQ,  at  kx=  0,  i.  e.  ,  where  the  dispersion  curve  crosses  the 
ky  axis.  It  will  first  be  shown  that  in  the  first  quadrant,  k^  is  a  single -valued,  mono- 
tonically  increasing  function  of  k  .  These  two  facts  indicate  that  the  inverse  function, 
ky=ky(kx),  is  single-valued  and  monotonically  increasing  in  the  first  quadrant  as  is 
depicted  in  Fig.  5.  Other  fundamental  properties  of  that  portion  of  the  dispersion  curve 
in  the  first  quadrant  of  Fig.  5  that  will  be  established  are:  1)  dk^/dk^  =0  at  kx=  0; 

2)  asymptotically  as  k^-*00,  k^k^/JX-l  and  the  dispersion  curve  everywhere  lies  above 
the  asymptote  k^  =  k^/  J X  -  1  ;  3)  the  value  of  k^  at  k^  =  0,  as  well  as  the  slope  of  the 
asymptote,  increase  with  uu .  One  question  that  has  not  yet  been  answered  analytically  is 

whether  the  surface  wave  dispersion  curves  have  inflection  points. 

To  see  that  in  the  first  quadrant  of  Fig.  5 ,  k^  is  a  single-valued  function  of  k^, 

observe  that  for  3>0,a>0  each  branch  of  -0  cot  8d  in  Fig.  3  and  each  branch  of  0tan0d 

2  2  2  2 

in  Fig.  4  intersects  the  circle  a  +  0  =  X(k^-kQ)  only  once.  Thus  for  a  given  and 

for  each  value  of  k  and  uu  there  will  be  only  one  set  of  values  (0,  a)  for  each  mode  and 

y 

hence  from  (23)  only  one  value  of  k^>  0  for  each  mode.  Therefore,  in  the  first  quadrant 

of  Fig.  5,kx  is  a  single-valued  function  of  k^.  That  kx  is  a  monotonically  increasing 

function  of  k  can  be  inferred  from  the  sign  of  dk  /dk  .  Since  k  and  k  satisfy  the 
y  °  x  y  x  Y 


3 


Fig.  5 


Dispersion  Curves  for  a  Typical  Surface  Wave  Mode  With  uu  as  a  Parameter. 
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surface  wave  dispersion  relation  D  (k  ,  uu)  =  0,  dk  /dk  for  fixed  uu  is  given  by 
i  s  t  x  y 


BD  /a  D 

dk  /dk  =  -  ~rT  /  --r— 

x  y  Bky/  akx 

defined  in  (11)  and  (14),  it  is  found  that 


Using  Dg  as  given  in  the  left-hand  side  of  (18),  with  a  and  8 


dk  k  ad(X-l)  +  k2/(k2-k2) 

_ x  _  V  x/  y  o' 

dk  "  TT  "  ad  +  1 

y  x 


(24) 


From  (24)  we  see  that  in  the  first  quadrant  of  Fig.  5,  dk^/dk  >  0  and  hence,  k^(k  )  is 
a  monotonic ally  increasing  function.  Furthermore,  (24)  shows  that  dk  /dk  =  0  at  k^=  0 
as  is  depicted  in  Fig.  5. 


As  k 


2  2  2  2 

the  value  of  8  at  the  intersection  of  the  circle  a  +  8  =  X(ky-kQ) 


and  any  one  branch  of  -8  cot  8d  in  Fig.  3  or  any  one  branch  of  8  tan  8d  in  Fig.  4 
approaches  a  constant.  Thus,  since  kQ  has  been  assumed  constant,  a  ~  Xky  as  ky-»® 
and  hence,  from  (23),  k^  in  the  first  quadrant  of  Fig.  5  is  asymptotically  given  by 
k  ~k  Jx-l  or  conversely  k  ~k  /i/X- 1 .  That  the  dispersion  curve  lies  above  the  asymp- 

x  y  . .  y  x 

tote  line  k^k^/VX-1,  as  shown  in  the  first  quadrant  of  Fig.  5,  can  be  deduced  from  the 
definition  of  8  given  in  (14).  Since  8  is  real  for  the.  surface  waves  and  X>  1, 

(X-l)k2-  k2  =  8^+k2(X-l)  >  0  and  therefore  in  the  first  quandrant  k  >  k^/  JX-1,  which 
proves  that  the  dispersion  curve  lies  above  the  asymptote  line.  When  uu  increases  but 
remains  below  uu^,  X  decreases  to  unity  and  hence  the  slope  of  the  asymptote,  1/JX-l, 
increases  as  is  depicted  in  Fig.  5.  Furthermore,  as  uu  increases  the  slope  of  the  line 
a  =  8/  Vx-l  in  Fig.  3  and  Fig.  4  increases.  Hence  the  values  ot  8  and  a  at  k^=  0,  as 
determined  from  the  intersection  of  the  line  a=  8/ Vx- 1  with  any  branch  of  -8  cot  8d  in 
Fig.  3  or  of  8tan8d  in  Fig.  4,  must  increase.  Because  k  increases  with  uu  and  X 


12  2  2 
decreases,  the  quantity  -^(a  +9  )  +  k 


at  k  =  0  increases  with 
x 


uu . 


k  must  increase  and  thus  the  magnitude  of  k 

°  y  y 

The  above-described  variation  with  uu  of  k  at  k  =  0  is 

y  x 


depicted  in  Fig.  5. 

Thus  the  fundamental  properties  previously  stated  for  the  surface  wave  dispersion 

curves  of  any  one  mode  are  seen  to  hold.  These  properties  indicate  that  the  dispersion 

curves  will  have  the  form  depicted  in  Fig.  5,  with  the  possible  exception  of  inflection 

points,  for  two  different  frequencies.  From  Fig.  3  and  Fig.  4  it  can  also  be  seen  that 

surface  waves  exist  for  all  ju  in  the  range  0  <  uu<  uu  .  Lastly,  since  in  the  first  quadrant 
dk  P 

of  Fig.  5 ,  >  0,  which  follows  from  (24),  and  since  the  dispersion  curve  for  uu*u>2  lies 


above  that  for  id  =  oij  '^2’  x  comPonent  of  ^  must  everywhere  be  negative.  That 

the  dispersion  curve  for  in  =  oj^  lies  above  that  for*  uj  =  tu^  <  ^3  ^°^ows  from  the  fact  that 

at  k^=  0  the  u)=uUj  curve  lies  above  the  uu—  uj^  curve  and  the  two  curves  never  cross 

since  V  u),  which  is  given  in  (42),  is  never  infinite.  The  observation  that  x  •  V,  uj>0  is 
-t  0  -t 

confirmed  by  the  analytic  expression  for  V  m  given  in  (42)  and  indicates  that  the  surface 

-t 

waves  are  of  the  backward  wave  type  with  respect  to  the  x  direction. 


GROUP  VELOCITY  AND  ENERGY  TRANSPORT  VELOCITY 

Having  established  the  basic  properties  of  the  dispersion  relation  of  the  surface 
waves  on  a  uniaxial  plasma  slab,  the  equality  of  group  velocity  and  energy  transport 
velocity  for  these  surface  waves  will  be  verified  by  direct  calculation.  This  relation  is 
given  in  Reference  (1)  as 


wO  1  00 

7kt‘"  =  J  £dz/  J 


where  s  represents  the  real  part  of  the  complex  Poynting  vector  ExH  and  w  the  time 
average  stored  energy  density.  To  this  end,  the  relative  field  amplitudes  are  first  cal¬ 
culated.  Since  one  of  the  coefficients  A^,  A^,  B1  and  is  arbitrary,  for  simplicity  let 

A1  =  -Aq(8-  j  a)  e^d  (26) 

where  Aq  is  arbitrary.  Then  from  (17)  it  is  found  that 

A2  =  -Ao(£!+  ja)e"j8d  (27) 


while  from  (16-a) 


B  =  -20A  e" 
1  o 


Using  the  dispersion  relation  in  the  form  given  in  (18),  which  is  valid  for  both  open-circuit 
and  short-circuit  bisection  modes,  it  follows  that 
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B3=  -26Ao  |41f  ='j28d«ad  ■  (29) 

With  these  expressions  in  equations  (12),  (13)  and  (15)  for  the  fields  in  the  three  regions, 
£  in  Region  1  is  found  to  be 

s=  4|A  |202u)€  e‘2a(z"d)[x  k  (k2-k2)  +  v  k  (a2+k2)]  (30) 

—  'o'  o  — o  x  o  y  -^-o  y  x 

while  in  Region  3  it  is 

£  *  4|Ao|Vweoe2“(z+d>  l^k^-k2)  (l,y.V)]  (3D 

and  finally  in  Region  2  it  is 

i  =  |Aol2«’Co(a2+B2l{2[iokx(kJ-^)+Ioky(02+^)]  + 


(32) 

00 

The  quantity  S  =  J  sd  z  is  now  calculated  to  be 
-  00 

S  =  4|Ao|2«Co(a2+3Z)|-xokx(k2-k2)[i+d] 

tioky[ik*  +  d(X-1,,ky-ko>]}  •  (33» 

In  order  to  determine  the  stored  energy,  observe  that  in  the  plasma  slab 


’  3  d  e  " 

3  in  e ' 

'1  0  0" 

=  € 

=  € 

0  1  +X  0 

3uu 

o 

3  ID 

O 

0  0  1 

(34) 


Thus  in  Region  2,  the  time  averaged  stored  energy  density, 


which  is  given  by 


(7,  8) 


is  found  to  be 


3  d  e 


J 


(35) 
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w 


=  |i*0i2^V)  - 


.cf±||ejP,3(z-d)+f^ 

(36) 

Outside  the  slab,  the  time  averaged  stored  energy  density  has  the  form 


*  ■  iU0l£l2+u„iH  I2] 


(37) 


so  that  in  Region  1 


w 


=  2 1 A  | 2  e  82[(k2+k2)(k2+a2)  +(k^-k2)2]  e‘2a(z'd) 
1  o'  o  o  y  x  '  y  o 


(38) 


while  in  Region  3 


w  =  2|Ao|2co82L(k2  +k2)(k2+a?)  +  (k‘-k‘)‘]  e‘ 


2  ,  2.2-,  2a(z+d) 


(39) 


Calculating  W  =  J  wd  z  ,  it  is  found  that 


W 


=  4 |A  1 2  e  (a2+e2)  ld(k2-k2)(Xk2-k2)  +  -i[k2k2  +(k2-k2)2]  } 
1  o1  o  l  y  o  y  o  a  x  y  y  o  J 


(40) 


In  deriving  the  above  power  and  energy  formulas,  extensive  use  has  been  made  of  the 
dispersion  relation  given  in  (18)  and  the  formulas  (11)  and  (14)  fcr  a  and  0.  Using  the 
above  expressions  for  S  and  W,  the  energy  transport  velocity  is  seen  to  be 


S  -x  k  (k2-k2)(- +  d)  +  y  k  [- k2+d(X-l)  (k2-k2)] 
—  _  — o  x'  y  o"a  "*"0  v  a  x  y  o  J 

W  -  * - 


y  °'  a  '  -o  y-a  x  •  •  •  y 

2  ,  2WV1  2  ,2,  ,  1  r  2,  2  ,  „  2  ,  2V2- 


d(k  -kl(Xk  -k‘)  +  -[k&k&  +  (k  -  k  )  ] 
y  o  y  o  a  x  y  y  o  J 


(41) 

3D 


In  order  to  compute  the  group  velocity  7  uj  ,  the  formula  V  n;  -■  -  7  D 

iS  t  iS  ^  ®  o  (JU 

from  implicit  function  theory  will  be  used  where  the  function  ^(k^  uu)  is  the  left-hand 


18 


side  of  (18).  It  is  found  that 


k2  X 


\  U)  =  U> 

-t 


i  rv  j\ 

•  x  k  (--  +  d)  +  y  k  [d(X-l)+ - ,  X  .  ] 

—  o  x  a  i-o  y  .  ?. ,  .2 ,  J 

a(g  +8  ) 


d(Xk2-k2)  +1  (k2-k2)  + 
y  o  a  y  o 


w 

_ 2LJL_ 

a(k2-k2) 
y  o' 


(42) 


2  2 

If  Loth  the  numerator  and  denominator  of  the  above  expression  are  multiplied  by  (k  -k  ) 

2  2  2  2  y  ° 

and  it  is  recognized  that  X(k^-kQ)  *  a  +3  ,  ^  tu  will  be  seen  to  be  identical  with  S^/W, 

as  predicted  in  Reference  1.  ~  * 

The  example  worked  out  above  also  illustrates  the  fact  that,  in  general,  the 

direction  of  £  as  well  as  its  magnitude  can  vary  with  z.  This  can  be  seen  from  equation 

2  Z  2  2 

(32)  for  s  n  Region  2  if  it  is  recognized  that  the  vectors  [x  k  (k  -  k  )  +  y  k  (8  +k  )] 

2  2  2  2  -o  x'  o  y'  y'K  x 

and  [xq  (kQ-  k^)  +  k^(k^-  8  )]  are  parallel  only  for  kx=0.  Since  the  coefficient  of 

the  first  vector  is  independent  of  z  while  the  coefficient  of  the  second  vector  depends  on 

z,  the  direction  of  the  vector  sum,  which  gives  £,  will  depend  on  z  for  all  surface  waves 

for  which  k  i-  0. 

x 
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APPENDIX 

The  purpose  of  this  Appendix  is  to  investigate  the  possible  contribution  to  a  sur¬ 
face  wave  field  frcm  the  plane  wave  fields  in  the  plasma  that  are  associated  with  the 
solutions  x  =  i^/k^-  of  the  plasma  plane  wave  dispersion  rela  ion.  The  vector  char¬ 
acter  of  these  plasma  plane  waves  is  that  of  H-type  modes  and  has  the  form 


(43) 


with  C  an  arbitrary  constant.  The  vector  character  of  those  plane  waves  in  the  air 
regions  that  have  the  H-type  mode  form,  and  will  thus  allow  a  simple  application  of  the 
continuity  conditions  at  z  =  ±d,  is 


e '  =  D[x  h  -  z  k  ] 

—a  — o  a  — o  x 


h'  =  — —  t-x  k  k  +  y  (k^-k^)-z  x  k  ] 
-a  WU-  -o  y  x  -^o'  o  y  -o  a  v 


(44) 


with  Ha  =  ±j  a  and  a  as  defined  in  (11).  Note  that  x  =  ±ja  also. 

In  Region  1,  xa  must  be  taken  as  -ja  to  ensure  that  E  and  H  are  zero  at  z  =  ®. 
Similarly,  in  Region  3,  xa  must  be  taken  as  ja-  Denoting  the  amplitudes  in  Regions  1  and 
3  as  Pj  and  D^,  respectively,  and  letting  Cj  and  be  the  amplitudes  of  the  plasma 
plane  waves  corresponding  to  x  =  -ja  and  x  =  ja  »  respectively,  the  continuity  conditions 
at  z  =d  result  in  the  equations 


-D1 e~ad  =-C1 e‘ad + C2 ead 
D: e"ad  =  CL  e'ad  +  C2  ead 


when  the  fields  in  Region  i.  are  assumed  to  be  the  sum  of  the  two  H-type  plane  waves.  The 
continuity  conditions  at  z  =  -d  can  be  written  as 


D3e'ad  =  -Ciead  +C2e'ad 
D3e-ad  =  Ciead+C2e-ad 


(46) 
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These  equations  have  only  the  trivial  solutions  C^=  C^=  D^=  D^=  0  and  hence  no  surface 
wave  can  exist  whose  fields  in  the  plasma  are  a  sum  of  the  two  H-type  plasma  plane  waves, 
which  propagate  as  x  =  ±ja  . 

The  physical  reason  why  no  surface  wave  exists  that  contains  the  above-mentioned 
plane  waves  is  that  the  waves  of  this  polarization  do  not  "see"  the  plasma,  since  the 
infinite  D.  C.  magnetic  field  along  y  prevents  the  electrons  from  moving  in  response  to 
an  R.  F.  electric  field  that,  as  in  this  case,  is  purely  transverst  to  y.  In  effect,  for 
waves  of  this  polarization,  no  slab  on  which  to  have  surface  waves  is  present. 
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